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We compute the absorption cross section and the total power carried by gravitons in the evapo¬ 
ration process of a higher-dimensional non-rotating black hole. These results are applied to a model 
of extra dimensions with standard model fields propagating on a brane. The emission of gravitons 
in the bulk is highly enhanced as the spacetime dimensionality increases. The implications for the 
detection of black holes in particle colliders and ultrahigh-energy cosmic ray air showers are briefly 
discussed. 
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I. INTRODUCTION 

Models with extra dimensions have emerged as the most successful candidates for a consistent unified theory of 
fundamental forces. The most recent formulation of superstring theory suggests that the five consistent string theories 
and 11-dimensional supergravity constitute special points in the moduli space of a more fundamental nonperturbative 
theory, called M -theory [l ]. In models of warped or large extra dimensions, standard model fields but the graviton 
are confined to a three-dimensional membrane II The observed hierarchy between the electroweak and the 
gravitational coupling constants is naturally explained by the largeness of the extra dimensions. 

The presence of extra dimensions changes drastically our understanding of high-energy physics and gravitational 
physics. Gravity in higher dimensions is much different than in four dimensions. For example, black holes with a 
fixed mass may have arbitrarily large angular momentum Q. The uniqueness theorem does not hold, allowing higher¬ 
dimensional black holes with non-spherical topology. A naive analysis suggests that higher-dimensional black holes 
should evaporate [B| more quickly than in four dimensions, due to the larger phase space. Moreover, brane emission 
should dominate over bulk emission because standard model fields carry a larger number of d.o.f. than the graviton 
I. However, a black hole does not radiate exactly as a black body and the emission spectrum depends crucially on 
the structure and dimensionality of the spacetime. A large graviton emissivity could reverse the above conclusion; if 
the probability of emitting spin-2 quanta is much higher than the probability of emitting lower spin quanta, the black 
hole could evaporate mainly in the bulk J7). A conclusive statement on this issue can only be reached if the relative 
emissivities of all fields (greybody factors) are known. This is particularly relevant for the phenomenology of black 
holes in models of low-energy scale gravity II Unequivocal detection of subatomic black holes in particle colliders 
fioj and ultrahigh-energy cosmic ray observatories 0 12 is only possible if a consistent fraction of the initial black 
hole mass [l3|, 13 is channeled into brane fields. 

If the center-of-mass energy of the event is sufficiently large compared to the Planck scale, quantum gravitational 
effects can be neglected and the black hole can be treated classically. It is commonly accepted that black holes with 
masses larger than few Planck masses satisfy this criterion. Under this assumption, if the fundamental gravitational 
scale is about a TeV, micro black holes produced at the Large Hadron Collider (center-of-mass energy = 14 TeV) and 
in ultrahigh-energy cosmic ray showers (center-of-mass energy > 50 TeV) can be considered classical. Throughout the 
paper we will assume that this is the case. (See Ref. 0 and references therein for a discussion of the uncertainties 
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deriving from this assumption.) 

The relative emissivities per degree of freedom (d.o.f.) of a classical four-dimensional non-rotating black hole are 
1, 0.37, 0.11 and 0.01 for spin-0, -1/2, -1 and -2, respectively M- In that case, the graviton power loss is negligible 
compared to the loss in other standard model channels. Since brane fields are constrained in four dimensions, the 
relative greybody factors for these fields approximately retain the above values in higher dimensions. The emission 
rates for the fields on the brane have been computed in Ref. 113] ■ The graviton emission is expected to be larger in 
higher dimensions due to the increase in the number of its helicity states. The authors have recently calculated the 
exact absorption cross section, power and emission rate for gravitons in generic H-dimensions E2- (See also Ref. 
0 .) The purpose of this paper is to discuss these results in more detail. The graviton emissivity is found to be highly 
enhanced as the spacetime dimensionality increases. Although this increase is not sufficient to lead to a domination 
of bulk emission over brane emission, at least for the standard model, a consistent fraction of the higher-dimensional 
black hole mass is lost in the bulk. 

The organization of this paper is as follows. In Section [H] we fix notations and briefly review the basics of gravita¬ 
tional perturbations in the higher-dimensional black hole geometry. In Section fTTTI we derive the field absorption cross 
sections from the absorption probabilities. Details of this derivation are included in the Appendix. The low-energy 
absorption probabilities and the cross sections for spin-0, -1 and -2 fields in generic dimensions are computed in 
Section Hvl This derivation is valid for scalars, vectors and gravitons in the bulk and generalizes previous results {l6i |. 
In Section 0 we prove that the high-energy behavior of the cross section is universal and reduces to the capture cross 
section for a point particle. Numerical results for the total power and the emission rates for all known particle species 
are obtained in Section IVII Finally, Section IyTTI contains our conclusions. 


II. EQUATIONS AND CONVENTIONS 

The formalism to handle gravitational perturbations of a higher-dimensional non-rotating black hole was developed 
by Kodama and Ishibashi |19j] (hereafter, KI). In this section we briefly review their main results. 


A. Number of degrees of freedom of gravitational waves 

In four dimensions, gravitational waves have two possible helicities [20| . corresponding to the number of spatial 
directions transverse to the propagation axis. In generic D dimensions, gravitational waves can be described by a 
symmetric traceless tensor of rank D — 2, corresponding to the CD representation of SO(D — 2). The number of 
helicities is 

r (D-2)(D-l) D(D ~ 3) , N 

M=~ -y--~1= ^ • (2-1) 

This representation is decomposed into tensor (T), vector (V), and scalar (S') perturbations on the sphere S D ~ 2 . 
These components correspond to symmetric traceless tensors, vectors and scalars, respectively. Tensors and vectors 
are divercenceless on S D ~ 2 . The number of d.o.f. are 

Mr = — 1^ - (D — 2), Af v = {D - 2) - 1, Ms = 1 . (2.2) 

The perturbations are further expanded in tensor, vector and spherical harmonics on S D ~ 2 . The total number of 
helicity states is obtained by considering the multiplicities of these components. A massless particle of spin J in four 
dimensions has helicities +J and —J, corresponding to the projections of the spin along the direction of motion. 
These states provide a representation of the little group SO(D — 2) = SO( 2), i.e. the group of spatial rotations 
preserving the particle direction of motion. In four dimensions all massless particles have two helicities since all 
non-singlet representations of SO( 2) have dimension two. The description of spin in D dimensions proceeds similarly 
to the four-dimensional case. For instance, in five dimensions there are three directions orthogonal to the direction of 
motion. The little group is SO( 3) and the helicities of the 5-dinrensional graviton are 2, 1 ,0, — 1 , —2. For a general 
discussion, see Ref. fill . 


B. Metric and master equations for gravitational perturbations 

The metric of the higher-dimensional non-rotating spherically symmetric black hole is 0, E2 


ds 2 = — fdt 2 + / 1 dr 2 + r 2 dVL 2 D _ 2 , 


( 2 . 3 ) 
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where 


/= 1 - 


th 

r D -3 


(2.4) 


The mass of the black hole is M = (D — 2)fi£>_2 rfr/(167 t(?), where £Id- 2 is the volume of the unit (D — 2)-dimensional 
sphere with line element and Q is Newton’s constant in D dimensions. Withouth loss of generality, we choose 

rn = 1. This amounts to rescaling the radial coordinate and u>. Hereafter, u> stands for urn- 

Let us consider weak perturbations of the geometry itO) by external fields. If the field is a scalar, the contribution 
to the total energy-momentum tensor is negligible. The evolution for the held is simply given by the Klein-Gordon 
equation in the fixed background. The evolution equations for the electromagnetic Held were derived by Crispino, 
Higuchi and Matsas and more recently by KI in the context of charged black hole perturbations. According to 
the KI formalism, the perturbations are divided into vector and scalar perturbations. The gravitational evolution 
equations were derived by KI E3 following earlier work by Regge and Wheeler J24j and Zerilli j2E| in four dimensions. 
The gravitational perturbations are divided in scalar, vector and tensor perturbations. Tensor perturbations exist 
only in D > 4. 

The evolution equation for all known Helds (scalar, electromagnetic and gravitational) can be reduced to the second 
order diHerential equation 


d 2 T 

dr 2 


+ (w 2 - V)tf = 0, 


(2.5) 


where r is a function of the tortoise coordinate r*, which is defined by dr/dr * = f(r). With the exception of 
gravitational scalar perturbations, the potential V in Eq. (12.51) can be written as 


K = f(r) 


1(1 +D- 3) , (D-2)(D-A) (l-p 2 )(£>-2) 2 


4r 2 


4 r £>-i 


( 2 . 6 ) 


The constant p depends on the field under consideration: 


P = 


0 for scalar and gravitational tensor perturbations, 

2 for gravitational vector perturbations, 

2/(D — 2) for electromagnetic vector perturbations, 

2 (D — 3)/(D — 2) for electromagnetic scalar perturbations. 


(2.7) 


For radiative modes, the angular quantum number l takes the integer values 

0,1... for scalar perturbations, 
l = { 1,2... for electromagnetic perturbations, 
2,3... for gravitational perturbations. 


The potential for gravitational scalar perturbations is 

V = f 


Q(r) 


16 r 2 H(r) 2 


where 


Q(r) = (D - 2) 4 (D - 1)V + (D - 2)(D - 1) 

x {4[2 (D - 2) 2 - 3 (D - 2) + 4 ]m + (D - 2)(D - 4)(£> - 6 )(D - l)}x 2 
- 12 (D - 2 ){(D - 6 )m + (D - 2)(D - 1 )(D - 4 )}mx + 16w 3 + AD(D - 2)m 2 , 

H(r) = m+^(D-2)(D-l)x, m = 1(1 + D - 3) - (D - 2), x = 


,D -3 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


and l > 2 has been assumed. In four dimensions, the equations for scalar and vector perturbations reduce to the 
Zerilli and Regge-Wheeler equations, respectively. In the limit l —► oo, Eq. E3 shows the universal behavior 


d 2 ^ . 2 J 2 nt 

77+ w -/-j ® = 0, 
art r z 


l 


( 2 . 12 ) 


The above equation guarantees that all perturbations behave identically in this limit. The effective potential for 
tensor perturbations is equal to the potential of a massless scalar field in the higher-dimensional Schwarzschild black 
hole background 0|. 
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Assuming a harmonic wave e lujt and ingoing waves near the horizon, we impose the boundary condition 

'F(r) -» e~ iujr * ~ (r - i)-<"/(^-3) ; r r+ . (2.13) 

At the asymptotic infinity, the wave behavior includes both ingoing and outgoing waves 

T(r) -► Te~ iur * + Re iur * , r* -> oo . (2.14) 

The physical picture is a wave with amplitude T scattering on the black hole. Part of this wave is reflected back with 
amplitude |i?| 2 and part is absorbed with probability |A| 2 = (|T| 2 — |i?| 2 )/|T| 2 . 


III. FROM ABSORPTION PROBABILITIES TO ABSORPTION CROSS SECTIONS 


The rate of absorbed particles for a plane wave of spin s and flux 4> s is 



(3.1) 


where a s is the absorption cross section. In Eq. El we sum over all final states and average over the initial states. 
The total cross section is obtained by summing the absorption coefficients for each single mode Z, A.f, weighted by 
the multiplicity factors. For a scalar field, the result is 1 


cr s=0 (D, oj) = C s=0 (D,cc)^^s | A?f 

l 


(3.2) 


where the normalization factor is 

C s = 0 (D,u) 


2tt\ D ~ 2 
uj ) 


Ad-2 


(47r)C D - 3 V 2 r[(T) - l)/2] 


(3.3) 


and the multiplicities of the scalar spherical harmonics are 


Nis 


{21 + D — 3)(Z + D — 4:)\ 

/!(£>-3)! 


(3.4) 


Equation iTOl follows from the decomposition of a plane wave in spherical waves. In four dimensions, Eq. (IQ) gives 
the well-known result Nis = 2i + 1. 

Spherical harmonics and partial wave expansion are different for vector and tensor fields. This leads to different 
expressions for the normalization and multiplicity factors. (This fact has been overlooked in the literature. See, for 
instance, Ref. J9|.) 

Consider a vector field V^ in the transverse gauge. The decomposition of a transverse plane wave in spherical waves 
gives the factor (see the Appendix D and Ref. 1231 1 


C s =\D,u;) 


C S=0 (D, w) 
D — 2 


(3.5) 


Following KI, we first decompose in the divergence of a scalar field plus a divergenceless vector. These components 
are then expanded in scalar and vector spherical harmonics. The multiplicities of the scalar harmonics are given by 
Eq. 13.411 . The multiplicities of the vector harmonics are 


The cross section is 


Niv 


1(1 +D- 3) (21 + D- 3 )(l + D- 5)! 
(I + 1)!(-D — 4)! 


ct5 =i = C s= i (A w) £ ^Nis I A?! 1 1 2 + N lv | Alt 
l 


(3.6) 


(3.7) 


1 See, for example, Ref. for the four-dimensional case and the Appendix B and Ref. nm for its D-dimensional generalization. 
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Finally, let us consider the graviton field. The decomposition of a plane wave in spherical waves yields the normalization 
factor (see the Appendix D for details) 


C a = 2 (D,u) = d J_ 3) C s =°(D,co). (3.8) 

The gravitational perturbation is first decomposed in scalar, vector and tensor components on the tangent space of 
S D 2 . These components are then expanded in scalar, vector and tensor harmonics. The corresponding multiplicities 
are Nis, Niy and 0] 


_ 1 (D- 1 )(D - 4 )(l + D- 2 )(l - 1)(2 1 + D- 3 )(l + D - 5)\ 
lT ~ 2 (Z + l)!(D-3)! 

respectively. The graviton absorption cross section is 

^=2 = C s=2 (D,w)^ [Nis IAf f 2 | 2 + iV* y | Afv 2 \ 2 + N lT \At ? 2 1 2 " 

l 


(3.9) 


(3.10) 


Note that in D = 4, accidentally, the conversion factors do not depend on the spin. The four-dimensional cross 
sections read 


a s (4,io) = ^y2(2l + l)\At\ 2 , s = 0,1,2. 


(3.11) 


(See also the discussion in Section fVl below .) 


IV. LOW-ENERGY ABSORPTION PROBABILITIES 


Since low frequencies (w<l) give a substantial contribution to the total power emission, it is worthwhile to derive 
an analytical expression for this limiting case. The low-energy approximation uses a matching procedure to find a 
solution on the whole spacetime for any value of p. Let us define the near-horizon region r — 1 -C 1/w and the 
asymptotic region r 1. To solve the wave equation in the near-horizon region, we set 

vfr = r^ 0 / 2 5 , (4.1) 

where /3° = {D — 2)(p + 1). With this substitution, the wave equation (12.511 is cast in the form 


r /3° 


30 dr (jr 0O d r s) 
Changing variables to 


2 , / ( ,n , n Q 3 , p(D-3)(D-2) , p 2 {D~2f 


u) H —~ ( —1(1 + D — 3) + 


+ : 


= 0. (4.2) 


v = 1 — 


1 


,D -3 ’ 


Eq. 14.21 reads 

(1 - v) 2 v 2 (3 - D) 2 d 2 E -(3- D)v( 1 - v) 2 ~ (3 - £>)) d v E + (( W r) 2 + va) \ 


where 


a = -1(1 + D - 3) + 


p(D -3)(D -2) p 2 (D — 2) 2 


2 4 

This equation can be put in a standard hypergeometric form by defining 

p(D — 2) u> a 


(4.3) 


= 0, (4.4) 


(4.5) 


Ph = 


D — 3 


D-3 


Uh = 


ah = 


(D - 3) 2 ’ 


(4.6) 
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and setting 


S = v Cl (l - v) C2 F, 

Cl = -iu h , 

1 . . lp°-l-(D- 3)6 

c 2 = -( 1 +^- 6 , = -- ^ 

1 

D~ 3 


b = J (1 + p h ) 2 - Aa h - 4u>l = ——-\/{D + 21 - 3 ) 2 - 4w 2 . 


The result is 


where 


t>(l — v)d 2 F + (7 — v(l + a + P)) d v T — aPF = 0 , 


7 = 1- 2iu> h , a = - (1 - p h - 2iu) h - b) , (3 = - (1 + p h - 2iu> h - b) . 

The most general solution of Eq. ED in the neighborhood of v = 0 is 

T = AF(a,P,"f, v) + B v 1-1 F{a — 7 + 1 ,/? — 7 + 1,2 — 7 , u). 


(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


Since the second term describes an outgoing wave near the horizon, we set B = 0. The asymptotic behavior of the 
near-horizon solution is 


U r -2ic-(0°-i)-(D-3)b/2 r(7)r(q + /?- 7 ) ( go_ 1)+(D _ 3)b/2 r(7)T( 7 -a-p) 

T(a)T(/3) + T( 7 — a)T( 7 — p) ’ 

where we have used the property of the hypergeometric functions: 

F(a,P, 7, v) = (l-v)"i - a -0 x r( 7 r[i)"r l (!g)~ 7) ~ ~ -a-P+ 1 , l-v) 

+ ^(^^^(7-/3) —7+a+/?+l, l-v). 

In the asymptotic region, Eq. iTQl) can be written as 


„ 2 _ 

o r zi i H- d r zi 

r 


The solution of this equation is 


j+ff-ni+D- 3 )+ ^-y- 2 ) + ^- 2); 


3 = C^ 1 -^ 2 J i/2 (ur) + C 2 r^y 2 Y l/2 (ur), 


= 0 . 


b/2 V 


where b = ^/(l — P 0 ) 2 — 4a. Expanding S for small wr, we obtain 

-c 2 


~ ~ c (^/ 2 ) b/2 r a-0°+b)/2 _ r (2/^) b/ 2 E(6/2) j 1 _ 0 o_~ b) /2 

1 T(1 + S/2) 


Equation 14.1811 is matched to Eq. 14.1411 . As ( D — 3)6 ~ b for u> <C 1, we find 


T( 7 - a - P)T{a)T(P) _ Ci n(%) b 


r(a + P- 7 )T (7 - a)T( 7 - P) C 2 r(l + |)T(|) ’ 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


Ci _ ( 2\ D+2l ~ 3 T(l + (l + ^=2) r ( 7 - a - /3)T(a)T(/3) 

C 2 ~ 


UJ 


7rT(a + P - 7 )r(y - a)r(y - P) 


(4.20) 


Cl = _ (2\ D+2l ~ 3 r (l + ^¥) 2 ( l + r (l + 2&) r (/3 ) 2 (l + 7&) sin ^ sin 
C 2 ~ l w 


21tt 
D -3 


7t 2 q: 2 sin7raT(—a) 2 


(4.21) 
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It is straightforward to show that the reflection coefficient is 


R Ci - iC 2 
T _ Ci + iC 2 ' 


(4.22) 


Therefore, the absorption probability in the u< 1 approximation is 


l-4 | 2 


1- \TZ\ 2 =4?r 


/WN D+21—2 

V2/ 


r 


C , 2l+p{D-2)\ 2 r C , 2l—p(D—2) \ 2 

2(D-3) J 1 ^ + 2(D-3) J 

r(i + i^) 2 r(/ + ^) 2 


(4.23) 


The result for the scalar waves M are obtained by setting p = 0 and l = 0,1.... Setting p = 0 (2) and Z = 2,3... 
we obtain the low-energy absorption probability for the gravitational tensor (vector) perturbations. The gravitational 
scalar perturbation cannot be dealt analytically. Numerical results give 


Pgrav scalar ~ 2 + 0.674L>- 0 ' 5445 . (4.24) 

Using Eqs. 63 , iron , dira . mu and (THiHl) it is straightforward to compute the low-energy absorption cross section 
for the gravitational waves. For instance, recalling that the four-dimensional absorption probabilities of scalar and 
vector perturbations are equal l30l and the tensor contribution vanishes, the contribution of the l = 2 mode in four 
dimensions is 


47T 4 

<n= 4S U 


This result agrees with the well-known result of Ref. 0. 


(4.25) 


V. HIGH-ENERGY ABSORPTION CROSS SECTIONS 

The wave equation can also be solved in the high-energy limit, where the absorption cross section is expected to 
approximate the cross section for particle capture. This has been verified in a number of papers for the scalar field. 
(See Ref. @ and references therein.) We will now extend this result to spin-1 and -2 fields in any dimension D. In 
the high-! limit, the multiplicities satisfy the relations 

Ni ,s — d _ - [Ni t s + N t >v ] = d _ ^ [IV; ,5 + Np V + Ni |T ] , Z —> oo , D> 2. (5.1) 

The above relations also hold for any l when D = 4. In that case, Chandrasekhar |30| showed that Ais = Ai v for 
any l and Nit{ 4,2) = 0. Therefore, in four dimensions the cross section does not depend on the particle spin. 

Since high frequencies can easily penetrate the gravitational potential barrier, the absorption probabilities of all 
fields, \A(D, l, u>)\ 2 , tend to 1 as u> —> oo. A WKB analysis shows that the absorption probability is nonzero for 
l < uj. Therefore, the cross section in the high-energy limit must include the contribution from all l < u>. The largest 
contribution to the cross section when to — > oo is given by high-! modes. Since the high-! limit of the wave equation 
is independent of the kind of perturbation (see Section nm the wave equation takes a universal form in this limit, 
i.e. Ais = Aiv = Ait- From Eqs. 13.51) - (13. 1 01) and (15.111 it follows that the absorption cross sections for spin- 0 , -1 
and -2 fields is equal in the limit u> —> oo. 


VI. TOTAL ENERGY EMISSION 


The total energy flux for gravitational waves is 


dE 

dt 


dEg dEy dEj' 

dt dt dt 



dio lu 
27 t e u> ! TH — 1 


(N^lAtfl 2 + N lv \A^ 2 \ 2 


Azt|A? t ~| 2 ) , 


( 6 . 1 ) 


where the Hawking temperature is Th = {D — 3)/(47t) and the counting of helicities is included in the multiplicity 
factors. The total absorption probabilities can be computed numerically. Equation (12.51) is integrated from a point 
near the horizon (typically r — 1 ~ 10~ 6 ), where the field behavior is given by Eq. (12.131) . The numerical result is then 













compared to Eq. (I2~T1 at large r. A better accuracy is achieved by considering the next-to-leading order correction 
terms jTL| 


where 


'F(r) 


T 


K) 


*(-;) 


r* 


oo. 


( 6 . 2 ) 


«l + D-3) + (D-2)jD-4)/A 

2u> v ’ 

The emission rates and the total integrated power for various fields are summarized in Tables nun For sake of 
comparison with previous works, the values for lower-spin fields are taken from Ref. jjjj]. (We checked these results with 
our numerical codes and found agreement within numerical uncertainties.) The results in the tables are normalized 
to the four-dimensional values. In four dimensions, the radiated power V is 

r° = 2.9xr 4 q 2 , 1P S=1/2 = 1.6 x 10 _4 r+ 2 , V s=1 = 6.7 x 10" 5 r+ 2 , V s = 2 = 1.5 x 1CT 5 r+ 2 . (6.4) 

The spin-0, -1/2 and -1 values are per d.o.f., whereas the graviton value includes the contribution of the two helicities 
2 . The four-dimensional emission rates are 


TZ S=0 = 1.4 x 10“ 3 r + 1 , K s = 1/2 = 4.8 x 10“ 4 , TZ S=1 = 1.5 x 1CT 4 r+ 4 , K s = 2 = 2.2 x 10~ 5 A/ 1 . (6.5) 

The above values for fermions, bosons and graviton agree with Page’s results (See Table I in Ref. |ljj). Some features 


TABLE I: Total radiated power V into different channels. The first three rows correspond to fields propagating on the brane. 
The last row gives the power radiated in bulk gravitons normalized to the four-dimensional case. 


D 

4 

5 

6 

7 

8 

9 

10 

11 

Scalars 

1 

8.94 

36 

99.8 

222 

429 

749 

1220 

Fermions 

1 

14.2 

59.5 

162 

352 

664 

1140 

1830 

Gauge Bosons 

1 

27.1 

144 

441 

1020 

2000 

3530 

5740 

Gravitons 

1 

103 

1036 

5121 

2 x 10 4 

7.1 x 10 4 

2.5 x 10 5 

8 x 10 5 


of the numerical results of Tableware worth discussing: 

• The relative contributions of the higher partial waves increase with D. For instance, in four dimensions the 
contribution of the l = 2 mode is two orders of magnitude larger than the contribution of the l = 3 mode. 
The l = 2 and l = 3 contributions are roughly equal in D = 8. More energy is channeled in l = 3 mode than 
in l = 2 mode for D > 9. (The largest tensor contribution in ten dimensions comes from the l = 4 mode.) 
Contributions from high l are needed to obtain accurate results for large D. For instance, in ten dimensions the 
first 10 modes must be considered for a meaningful result. Therefore, precise values for very large D require the 
most CPU-time. The values in Table [I] have a 5% accuracy. 

• The total power radiated in gravitons increases more rapidly with D than the power radiated in lower-spin 
fields. This is due to the increase in the multiplicity of the tensor perturbations, which is larger than the scalar 
multiplicity by a factor D 2 at high D. Therefore, the main contribution to the total power comes from tensor 
(and vector) modes. For instance, in ten dimensions the tensor modes contribute roughly half of the total power 
output. 

Table HT1 gives the fraction of radiated power per d.o.f. normalized to the scalar field. In four dimensions, the graviton 
channel is only about 5% of the scalar channel. Therefore, the power loss in gravitons is negligible compared to the 
power loss in lower-spin fields. This conclusion is reversed in higher dimensions. For instance, the graviton loss is 
about 35 times higher than the scalar loss in D = 11. Graviton emission is expected to dominate the black hole 
evaporation at very high D. The particle emission rates per d.o.f. are shown in Table III. The relative emission 


2 Since the number of graviton d.o.f. (helicities) depends on the spacetime dimensionality D, here and throughout the paper we give the 
total rate and the total power for the graviton field, rather than the rate and power per d.o.f. 
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TABLE II: Fraction of radiated power per d.o.f. normalized to the scalar field. The graviton d.o.f. (number of helicity states) 
are included in the results. 


D 

4 

5 

6 

7 

8 

9 

10 

11 

Scalars 

1 

1 

1 

1 

1 

1 

1 

1 

Fermions 

0.55 

0.87 

0.91 

0.89 

0.87 

0.85 

0.84 

0.82 

Gauge Bosons 

0.23 

0.69 

0.91 

1.0 

1.04 

1.06 

1.06 

1.07 

Gravitons 

0.053 

0.61 

1.5 

2.7 

4.8 

oo 

00 

17.7 

34.7 


TABLE III: Fraction of emission rates per d.o.f. normalized to the scalar field. The graviton result includes all the helicity 
states and counts as one d.o.f. 

D 

4 

5 

6 

7 

8 

9 

10 

11 

Scalars 

1 

1 

1 

1 

1 

1 

1 

1 

Fermions 

0.37 

0.7 

0.77 

0.78 

0.76 

0.74 

0.73 

0.71 

Gauge Bosons 

0.11 

0.45 

0.69 

0.83 

0.91 

0.96 

0.99 

1.01 

Gravitons 

0.02 

0.2 

0.6 

0.91 

1.9 

2.5 

5.1 

7.6 


rates of different fields can be obtained by summing on the brane d.o.f. For instance, the relative emission rates of 
standard model charged leptons (12 d.o.f.) and a 11-dimensional bulk graviton are roughly 1:1. This ratio becomes ~ 
40:1 in five dimensions. To illustrate the relevance of these results for black hole in particle colliders, let us consider 
the minimal U(l) x SU(2) x SU{3) standard model with three families and one Higgs field on a thin brane with 
fundamental Planck scale = 1 TeV. For black holes with mass ~ few TeV the Hawking temperature is generally above 
100 GeV. The temperature of a six-dimensional black hole with mass equal to 5 (100) TeV is ~ 282 (133) GeV. The 
temperature increases with the spacetime dimension at fixed mass. Therefore, all d.o.f. can be considered massless. 
(Considering massive gauge bosons does not affect the conclusions significantly.) The spin-0, -1/2 and -1 d.o.f. on the 
brane are 4 (complex Higgs doublet), 90 (quarks + charged leptons + neutrinos) and 24 (massless gauge bosons), 
respectively. The relative emissivities for this model are shown in Tabic Hvl Although the graviton emission is highly 
enhanced, the large number of brane d.o.f. implies that the brane channel dominates on the bulk channel. However, 
power loss in the bulk is significant and cannot be neglected at high D; about 1/4 of the initial black hole mass is lost 
in the 11-dimensional bulk. This implies a larger-than-expected missing energy in particle colliders. 


TABLE IV: Percentage of power going into each field species for the minimal 1/(1) x SU(2) x SU( 3) standard model with three 
families and one Higgs field above the spontaneous symmetry breaking scale. The four-dimensional results are taken from Ref. 
fl5j and the higher dimensional results for fermions and gauge fields are taken from Ref. [|f. 


D 

4 

5 

6 

7 

8 

9 

10 

11 

Scalars 

00 

CD 

4.0 

3.7 

3.6 

3.6 

3.5 

3.3 

2.9 

Fermions 

83.8 

78.7 

75.0 

72.3 

69.9 

66.6 

61.6 

53.4 

Gauge Bosons 

9.3 

16.7 

20.0 

21.7 

22.3 

22.2 

20.7 

18.6 

Gravitons 

0.1 

0.6 

1.3 

2.4 

4.2 

7.7 

14.4 

25.1 


VII. CONCLUSIONS 

In this paper we have computed the absorption cross section and the total power carried by gravitons in the 
evaporation process of a higher-dimensional non-rotating black hole. We find that the power loss in the graviton 
channel is highly enhanced in higher-dimensional spacetimes. This has important consequences for the detection of 
microscopic black hole formation in particle colliders and ultrahigh-energy cosmic ray observatories, where a larger 
bulk emission implies larger missing energies and lower multiplicity in the visible channels. Despite the increase in 
graviton emissivity, for 4 < D < 11 a non-rotating black hole in the Schwarzschild phase will emit mostly on the 
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brane due to the higher number of brane d.o.f. 3 However, black hole energy loss in the bulk cannot be neglected in 
presence of extra dimensions. Graviton emission is expected to dominate the black hole evaporation at very high D. 
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APPENDIX: DERIVATION OF THE CROSS SECTIONS 

In this appendix we derive in detail the cross sections for spin-0, -1 and -2 fields in generic dimensions. The harmonic 
expansion is performed on the spherical manifold S n with metric 7 a b and coordinates y a (a = 1,..., n = D — 2). We 
also define the set of coordinates x l = (r, y a ). In the limit r —> 00 the wave vector of a massless particle is 

= (w, k), (A.l) 

where k = ( k 1 ) = tok. The direction of the plane wave is given by the (n + l)-dimensional unit vector k. 


Scalar perturbations 


Let us consider an incoming scalar plane wave with unit flux and wave vector k = ivk at infinity: 

-> . (A.2) 

This wave can be expanded in spherical waves \E , ( A,a ’) parametrized by a discrete index A: 

<#’“) = a Ck, w; A)tf (A - w) . (A.3) 


In general, A = (km). However, we leave the index implicit to simplify the generalization to higher spins. The 
asymptotic behavior of the fields is 


v[r(V“) £_ yW 

r n/2 1 5 


where the spherical harmonics, Y^\ satify the normalization and orthogonality conditions 


J (iH„V (A) V (A,) = 6 XX '. 


(A.4) 


(A.5) 


From the above equations it follows 

r n J dn n 'I>* (A ’ a ’ ) <F (A ’ w) = 1, (A.6) 

i.e. the wave 'ld A,aJ ) describes one ingoing particle per unit time. If A s=0 (A,w) is the absorption coefficient for a 
spherical wave with angular profile V^ A ^, the cross section associated to the plane wave ( k,u >) is 

= £ | a(k,u>; A)| 2 |A S=C V, A)| 2 . (A.7) 

A 


3 This conclusion holds only if the particle content at TeV energies is mostly made of fields propagating on the brane. If non-standard 
model d.o.f. open up at the TeV scale, such as vector multiplets propagating in the bulk, black hole emission could occur mostly in the 
bulk. 
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The average on the wave direction of the initial state is obtained by integrating k over the sphere and dividing by the 
unit volume Q n : 


= — Y 

o 


a(fc,o;;A)| 2 |A s =>,A)| 2 . 


From the previous equation, we obtain 


C s =°(n + 2,Lo) = —j dn^\a(k,cu;X)\ 
The coefficients a are given by the coefficients of the expansion 

0 ( *’ w) = ^a(fc,w; A)^ (a ’ w) , 


(A.10) 


where 


VJ>(A,oj) = e ^ y(A) = g-iwt-ifea 

are defined in the flat (n + l)-dimensional Euclidean space. From 


J ^+ 1 cc «r(A',u/) v j / (A,u;) = 27rJ ^ _ w ') ( 5 aa , 


it follows 


Finally, the normalization factor is 


C s =°(n + 2,w) = 


a(k,u> ; A)| 2 = 


1 / 2 tt 


V ^ 


(A.ll) 


(A.12) 


(2Z) 2 <5(w - w') J(w' - a/")|a(fc, w; A) | 2 = j d n+1 xd n+1 x'^ x ’^ ( ; x )<p’"') (s') * (A, “" ) (o'). (A. 13) 

Multiplying by u> ,n and integrating in dfl^dw', we obtain 

u> n (2Tr) 2 6{uj -w") j dn^\a(k,uj-,X)\ 2 

= J d" +1 :rd" + V#* (A ’ w) (£)^ (A ’ w ''V) J dw'u ,n J d^e^-^ 


(A.14) 


(A.15) 


Vector perturbations 


Let us consider an incoming transverse vector plane wave with unit flux and wave vector at infinity: 


p,fc,u;) ^ iut— ikx 


m / ’ 


(A.16) 


where (/c) are the transverse polarization vectors (p = 1 ,,n) and 


k l ef = 0 , efe^ 6 lj = S pp , ^ = dy + terms in ki, kj . 


(A.17) 
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We assume x 1 to be asymptotically Euclidean coordinates for simplicity. The wave can be expanded in spherical 
waves vj/( A >“) parametrized by a discrete index A, which includes the polarizations: 




(A.18) 


The asymptotic fields 4/^ are 


\T/(A,w) 


-iiot—icor 


yTlj 2 


\ rY 0 


(A) 


(A.19) 


where Ya^ satisfy the normalization and orthogonality conditions 

J dS2„Y a ^Y b (y) 7 oh = 6 AA ' . 

From the above equation it follows 

r n j dn n $*( A .")$( i *.‘-y*' = l. 


(A.20) 


(A.21) 


i.e. a wave V F ; ', A ’^ describes one ingoing particle per unit time. If A S_1 (A, u>) is the absorption coefficient for a spherical 
wave with angular profile Ya X \ the cross section associated to the plane wave (p, k,u>) is 


= E l“(P. w; A)| 2 |A S=1 ( W , A)| 2 . 


(A.22) 


The average on the wave direction and polarization of the initial state is obtained by summing on p, integrating k 
over the sphere and dividing by the unit volume and the number of polarizations n: 


1 


a = 


Tlilr 


E / ^«EK^;A)| 2 |^’ A) | 

A J p=1 


(A.23) 


From the previous equation, we obtain 


C s = 1 (n + 2,u) = 


71 fL 


4^I4,M;A)I 2 . 

P=1 


The coefficients a are the coefficients of the expansion 


where 


§( A ’ W ) = £ 


—icot—icor 


o 


= e 




(A.24) 


(A.25) 


(A.26) 


r n/2 ^ r Y a ( A) J ’ 

are defined in the flat (n + l)-dimensional Euclidean space. From 

J d n+ 1 x^* (A, ’ a ’ ,) ^‘ (A ’ a, ' ) “ 

it follows 

(2it) 2 5(uj - lu')5(u' -u")\a{p,k,u-,X)\ 2 = [ . (A.28) 


= 27T(f(a; — o/)c>aa' 


(A.27) 
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Multiplying by u ' n , summing over p and integrating in dCt^dw', we obtain 

r n 

u) n (2ir) 2 6(u> - u>") / df^ fc) ^|a(p,fc,w;A)| 2 


= J d n+1 xd n+1 x'^* {x ’ w \x)^^ u "\x') J dw’J n J dSifi f^e i(p) (k)e* (p) (k)j e -ft*-*) 

= (2 ir) n+1 J d n+1 xd n+1 x'6{:x - x')*? X ’“\x)& {X ’ w "\x') (<jj + terms in k\ kj) = (2Tr) n+2 S(uj - to"). 


Finally, the normalization factor is 


c '"‘ ( " + 2 -" ) = Jn f dn i. t| El«(p.*.“0)i 2 = Jr (§ 


(A.29) 


(A.30) 


Gravitational perturbations 

In our normalizations, a spin-two field has number flux (This choice is possible because the waves are 

monocromatic, i.e. ~ e IUJt .) If is multiplied by a suitable factor, it can be interpreted as a metric perturbation 
h ltv . The multiplication factor depends on to and h, but its explicit form is not relevant in the derivation of the cross 

section. An incoming transverse spin-two plane wave (f >)with unit flux and wave vector at infinity k = u>k is 


±(p,k,cj) _ —iujt—ikx 


0 0 \ 

o eUk) ) ’ 


where e p Ak) are the transverse traceless polarization vectors, p = 1 ,... AT = (n — l)(n + 2)/2 and 


(A.31) 


k% = 0 , (Zjf&S *# 1 = S pp ' , e iJ ePkl = 6 \A) ~ “M** + terms in k u kj,k k , k l . (A.32) 


This wave can be expanded in spherical waves 




The asymptotic fields are 


vu( A » a; ) 

ins 


o 

o 

0 

0 

to 

e"! 


where Y^ 1 satisfy the normalization and orthogonality conditions 

Jdn n Y^Y}pr c i bd = s xx ' 


From the above equation it follows 


r n J dn^Yu)^,*) nu = 1} 


(A.33) 


(A.34) 


(A.35) 


(A.36) 


i.e. a wave describes one ingoing particle per unit time. If A S_2 (A, u>) is the absorption coefficient for a spherical 

wave with angular profile Y^\ the cross section associated to the plane wave ( p , k,co) is 

a Pl,u _ ^ \a(p, k,w; A)| 2 |A s=2 (w, A)| 2 . 


(A.37) 
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The average on the wave direction and polarization of the initial state is obtained by summing on p, integrating k on 
the sphere and dividing by the unit volume Q n and the number of polarizations TV": 


M 


(7 = 


TVTL 


H Wp,^w;A)| 2 3 |ai (i 

A J p =1 


a;,A) |2 


(A.38) 


From the previous equation, we obtain 


C s ~ 2 (n + 2, ui) 


1 

7VT2„ 


W 

P=1 


(A.39) 


The coefficients a are the coefficients of the expansion 


where 




= /)°t(p,k, 


A 




A)^’“ } : 


’F 




g— iujt—iujr 

( 0 

0 

y>n/2 

v° 

PUX 


^rXu,) = e -i^t-ife e P (fc) 


(A.40) 


(A.41) 


are defined in the flat (?z + l)-dimensional Euclidean space. From 

J d n+1 x^ X '’ u, ' ) ¥ j ^ = 27t5(w - o/)<S AV , 


(A.42) 


it follows 

(2tt) 2 8(u — uj')8(lu' — u/^lafp, k,u>; A )| 2 


d n+1 xd n+1 x'^ X ’ u \x)^^’ u '\x)^* k <*' k ' u ' ) (x')V kllx ’ u "Xx'). (A.43) 


Multiplying by ui' n , summing over p and integrating in dCt^ dui ', we obtain 


M 


w"(27T ) 2 5{u-u") / |a(p, A:,w; A )| 2 

J P =l 


= / d n+1 a;cr i+1 x , ^*.^’ a,) (a:)4 ffci *' A, “"' ) (a; 7 ) 


dcu'u' n J dflfi \ J^e ij(p) W4 ( i P) W J e 

= (27T)- 1 / d n+ 1 xd n+ 1 a; , (5(f-f , )^ (A ’“ ) (a:)'F i;^ ( A '“ ,, )( a ; , ) + terms in k\k j ,k k ,k, 

= { 2 n) n+ 2 S{to - w "). 


Finally, the normalization factor is 


a r 


C s=2 (n + 2,w) = 


TVfi* 


A ) I 2 = 

p=i 


1 /2 tt\ ? 


AfCl n V U> J 


(A.44) 


(A.45) 
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